New Global Minima for Thomson's Problem of Charges on a 

Sphere 

Eric Lewin Altschuler 

Mt. Sinai School of Medicine 1425 Madison Avenue, 
Box 1240 New York, NY 10029 
email: Eric.Altschuler@mssm.edu 

Antonio Perez-Garrido 
Departamento de Fisica Aplicada Universidad Politecnica de Cartagena, 
Campus Muralla del Mar Cartagena, 
30202 Murcia, Spain 
email: Antonio. Perez@upct.es 

Abstract 

Using numerical arguments we find that for N = 306 a tetrahedral configuration (T/^) and for 
N = 542 a dihedral configuration (D5) are likely the global energy minimum for Thomson's problem 
of minimizing the energy of N unit charges on the surface of a unit conducting sphere. These 
would be the largest N by far, outside of the icosadeltahedral series, for which a global minimum 
for Thomson's problem is known. We also note that the current theoretical understanding of 
Thomson's problem does not rule out a symmetric configuration as the global minima for N = 306 
and 542. We explicitly find that analogues of the tetrahedral and dihedral configurations for N 
larger than 306 and 542, respectively, are not global minima, thus helping to confirm the theory 
of Dodgson and Moore (Phys. Rev. B 55, 3816 (1997)) that as grows dislocation defects 
can lower the lattice strain of symmetric configurations and concomitantly the energy. As well, 
making explicit previous work by ourselves and others, for N < 1000 we give a full accounting 
of icosadeltahedral configuration which are not global minima and those which appear to be, and 
discuss how this listing and our results for the tetahedral and dihedral configurations may be used 
to refine theoretical understanding of Thomson's problem. 



What configuration of unit cliarges on tlie (surface) of a unit conducting sphere min- 



imizes the Coulombic energy Y^i^jl/vij Beyond physics this question has utility in 
assembly of biological [2] and chemical structures 0, Q|, to benchmark optimization al- 
gorithms, and, as well, in mathematics Smale jsj has noted this question to be a Hilbert 
problem for the Twenty-First Century. For 2 < N < 100, the original question asked by J.J. 
Thomson a century ago, and a bit beyond, there is agreement of all numerical 



and theoretical ll| methods suggesting that the global minimum configuration has been 
found. However, for larger N owing to an exponential growth in good local minima 
finding global minima, general principles or insights for minimization, or even methods or 
cases to test hypotheses has proven extremely difficult. For = 10(/i^ + k"^ + hk) + 2 
highly symmetric icosadeltahedral configurations can be constructed. While it was initially 



br large N adding defects 



17| . Here we note a tetra- 



thought that such configurations might be global minima 
to the icosadeltahedral lattice lowers the energy 3, M , lli 
hedral configuration for = 306 and a dihedral configuration for = 542 which based on 
numerical arguments appear to be global minima, the largest such A^ by far, aside from the 
icosadeltahedral series, for which a global energy minimum configuration is known. Study of 
tetrahedral analogues larger than 306 and dihedral analogues larger than 542 help confirm 



the theory 18| on defects lowering lattice strain and energy. As well, we note that the fact 
that lattice configurations fail to be global minima for A^ > 800, may help explain why M0154 
anions self-assemble into a spherical superstructure with a non-lattice number of subunits 



, the presumed global minimum 
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111 ] . We had previously 



For A^ = 78, as originally suggested by Edmunson |11| 
configuration has tetrahedral (T/j) symmetry (Figure la)!?,! 
suggested that an analogue of this configuration with 306 charges (Figure 16, see below for 
method of construction of this analogue) also might be a global minimum and done limited 
numerical testing of this idea [12!]. Now, here we have extensively checked numerically on 
over one thousand runs, and have found no configuration of lower energy. (Numerically to 
look for non-lattice configurations with energies lower than the lattice configuration we used 
random initial configurations followed by a local steepest descent method. For A^ = 306 we 
could not find a configuration in one thousand runs with an energy lower than the lattice 
energy, though for many other A^ considered in this paper including A^ = 1218 and 4866, and 
smaller A^ with icosadeltahedral lattices (see Tables I and II, and Figure 2) configurations 



2 



with energies lower than the lattice energy could be found in '^^^ fifteen runs.) But one 
thousand runs is orders of magnitude less than the estimated [T] over 1.5 million local 
minima for = 306. As well, without an analytic proof we could not be certain the 
tetrahedral configuration in Figure lb is the one of minimum energy for A^ = 306. Current 
theoretical understanding of Thomson's problem jl3, 16, 17, 18, 1^, to be discussed below, 
does not exclude the possibility of a symmetric configuration as the global energy minimum 
for A^ = 306. Of course, our configuration stands open to challenge. 

Using a genetic algorithm Morris, Deaven and Ho confirmed previously found pre- 
sumed global minima for A^ < 112 and gave their likely global minima for A^ < 200. Fitting 
the energy of these minima for A^ < 200 to the function: E = (l + aN^^/"^ + bN~^^'^^, 

(see 13] and refs. therein for an explanation of why this function was used) for A^ >100 they 
found particularly deep minima with respect to this function for icosadeltahedral configura- 
tions for A^ = 122, 132 192 and for dihedral configurations for A^ = 137, 182 and 187 and 
D2 configuration for N = 146. By our numerical testing as well these dihedral configurations 
appear to be global minima, though we have no analytic proof and numerically we are orders 
of magnitude short [7] to even sample a majority of local minima. We next looked at the 
higher split analogues of these dihedral configurations with AN — 6 charges. (See below for 
the method of construction of such analogues.) The next larger dihedral analogues for A^ = 
146, 182 and 187-A^ = 578, 722, and 742 respectively-are found not to be global minima 
after only a few runs (E.L.A. and A.P.G. data not shown). However, the next larger dihedral 
(Ds) analogue of A^ = 137-A^ = 542-after over one thousand runs appears possibly to be 
a global minima. As for the tetrahedral configuration for A^ = 306, we have no analytic 
proof of this proposition and given the huge number of local minima for an A^ this large (71] 
our numerical runs only begin to address the question. Conversely, we easily found that the 
next analogue of A^ = 137, A^ = 2162, is explicitly not a global energy minimum. See Figure 
3 for dihedral configurations for A^ = 137, 542 and 2162. 

Euler's theorem asserts that when a (convex) polyhedron is constructed by joining points 
on a sphere the number of faces (F) plus the number of vertices (V^) is equal to the number 
of edges (E) + 2; F + V = E + 2. For A^ > 12 this has the result that in addition to the 
sixfold coordinated points of a planar two-dimensional lattice (hexamers), there must be 
at least twelve points of fivefold coordination (pentamers). The tetrahedral configurations 
for A^ = 78 and 306, dihedral configurations for A^ = 137 and 542, and icosadeltahedral 
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configurations for = 10(/i^ + k"^ + hk) + 2 liave exactly twelve pentamers and the rest 
hexamers (see, e.g., Figure la, b, c, e. Figure 2a, b and Figure 3a, b, c). Larger analogues 
{split configurations) of the tetrahedral configuration for N = 78 (for =306, 1218, 4886 
and for the dihedral configuration for = 137 (A^ = 542, 2162) (Figure 16, c, e and Figure 
36, c) are made as follows: in addition to the A^ charges preexisting place one charge at the 
center of each of the 3A^ — 6 edges. (If all the charges were sixfold coordinated hexamers 
there would be 3A^ edges, six must be subtracted from this to take into account the twelve 
fivefold coordinated pentamers.) Then relax to the final position by a local gradient method. 
The resulting configuration has A^ + (3A^ — 6) = 4A^ — 6 charges. (See Figure 16, c, e). Some 
icosadeltahedral configurations can be made (split) as analogues of smaller ones (see e.g. 
A^ = 1242 in Figure 26, a larger analogue of A^ = 312 in Figure 2a). A method for making 



121. 



icosadeltahedral configurations de novo has been discussed previously 

Conversely, to the cases of A^ = 78 and 306, and A^ = 137 and 542 for the larger analogues 
we have studied (A^ = 1218, 4886, 2162) we find that adding dislocation defects to the lat- 
tice produces a configuration with lower energy (Fig. Ic — / and Fig. 3c, d ). Similarly, for 
A^ > 792 icosadeltahedral configurations with dislocation defects, additional fivefold coordi- 
nated points and then necessary also sevenfold coordinated points (septamers), have lower 
energy than symmetric lattice configurations, while for smaller A^ the symmetric lattice con- 
figurations appear to be global minima (see Tables I and II, see Figure 2 for an example), 
though further numerical testing may show that some such configurations are not global 
minima. 

These numerical results on the A^ at which tetrahedral, dihedral and icosadeltahedral 
configurations fail to be global energy minima are in remarkable concordance with a theory 
given by Dodgson and Moore originally for icosadeltahedral configurations: Using con- 
tinuum elasticity theory [3| they studied the energy cost of a pair of pentamers, compared 
with a pure hexagonal lattice, and suggested that dislocation defects-extra fivefold coor- 
dinated points, with (necessarily) paired sevenfold coordinated points-would lower lattice 
strain and energy for A^ in the 500-1000 range. Similar reasoning should apply to the 
tetrahedral configurations in the N = 78 series and dihedral configurations in the A^ = 137 
series. Our results given here are strong confirmation of Dodgson and Moores theory 1^ . 

The fact that apparently for A^ > 792 all symmetric tetrahedral, dihedral or icosadeltahe- 
dral lattices are not global minima, along with the exponential growth in good local minima 
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may help explain why the number of M0154 anions which self-assemble into a spherical su- 

n 

perstructure |4] is a non-lattice number 1165, rather than, for example, 1172 the closest 
icosadeltahedral lattice or the tetrahedral lattice at 1218, while for small N self-assembly 
often produces a symmetric lattice configuration 

as. 

For icosadeltahedral configurations for < 792 whether or not a lattice configuration is 
a potential global minimum depends not only on the magnitude of A^, but also apparently 
on the details of the lattice itself. All lattice A^ are listed in Tables I and 11. As can be seen 
for A^ = 42, 92, 162, 252, 362, 432, 492, 572, 642, 732 the icosahedral lattice configuration is 
manifestly not the global minimum, while for the other A^, the lattice configuration appears 
to be so. 

In an icosadeltahedral lattice A^ = 10(/;.^ + k"^ + hk) + 2, to go from the center of one 
pentamer to the center of an adjacent pentamer one moves h steps along one basis vector 

o the lattice, and then k steps in the other. We noted previously [12| (also discussed in ref. 

3|) that the energy in an icosadeltahedral lattice configuration with a large ratio of h to 
k {h > k) may be increased due to the vertices of the pentamers being closely aliged (or 



lat as N grows the 



13|, 
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perfectly aligned in a lattice with k = 0). It has previously been noted t 
icosadeltahedral lattice configuration may not be the global minimum 
though we have not seen any explicit published accounting of the A^ for which the lattice 
fails to be a global minimum. This is given in Tables I and II (in addition to the rule 
that for A^ > 792 the lattice is not the global minimum). A clear pattern emerges: for an 
icosadeltahedral configuration with k = 0, besides the extremely exceptional case of A^ =12, 
the icosadeltahedral lattice configuration is not the global minimum. For the three small A^ = 
42, 92 and 162 the apparent global minimum configuration has exactly twelve pentamers, and 
no dislocation defects, but arranged in an non-icosadeltahedral configuration, likely lowering 
the energy cost of having the vertices of pentamers aligned. For A^ larger than 162 in the k 
= series, the apparent global minimum incorporates dislocation defects. For A^ = 432 and 
higher for lattices with k = 1, the lattice is also not the global minimum in accordance with 
the notion stated above that the energy of the lattice is increased by relative alignment of the 



vertices of the pentamers. The data in Tables I and II ma^. 
predictions for global enery minima: The current theory 



DC useful in refining theoretical 
liQco„ectlyp.ed,ct.thatfo. 



A^ > 1000 icosadeltahedral and tetrahedral configurations will not be global minima, but 
does not yet account for various cases for A^ < 1000. As well, for k = other theoretical 
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work 



predicts dislocations lowering the energy for > 300, but the first instance of 



this is found for = 252. 

For tetrahedral or dihedral lattices we have not yet been able to find an obvious rule or 
principle to predict for which the lattice configuration is a global minimum. Indeed, while 
all groups using a variety of different methods find the configuration a global minimum 
for N = 78 10], and we have found similarly for = 306, there are a number of A^ < 100 
for which Th configurations are not the global minima [10|, lul, and also we have found 
that the next analogue of the global minimum for A^ = 100 which has T (though not Th) 
symmetry, A^ = 394, is not a global minimum (E.L.A. and A.P.G., data not shown). And 
we see no obvious difference between the N = 78 lattice and the others to explain why not 
only for A^ = 78, but the next higher analogue the lattice appears to be the global minimum. 
Similarly, we have appreciated no obvious reason why the next D2 analogue of 146-578-or 
the next -D5 analogues of 182 and 187-722 and 742-are not global minima, but the next 



D5 ana. 
Moore 



ogue of 137-542-appears to also be a global minimum. The theory of Dodgson and 
13l ] does not predict a priori that defects would lower a lattice energy for an A^ of 306 
and 542, as it does for the next analogues 1218 and 2162. Perhaps, whatever still unknown 
reasons that explains the good minimum for the lattice for N = 78 and lattice for 
A^ = 137 also permit the A^ = 306 and 542 analogues to be global minima. 

We thank the anonymous reviewers for comments extremely helpful in revising the 
manuscript and also inspirational leading to finding a new possible global minimum. A.P.G. 
would like to aknowledge financial support from Spanish MCyT under grant No. MAT2003- 
04887. 
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FIG. 1: Tetrahedral configurations. For a) N = 78 {E = 2662.04647), b) N = 306 {E = 
43862.56978), the symmetric tetrahedral lattice configuration appears to be the global energy min- 
imum, while for N = 1218 c) {E = 718284.03747), d) {E = 718281.63110) and N = 4886 e) {E = 
11651484.51295), /) {E = 11651440.24177) addition of defects to the lattice produces a configura- 
tion of lower energy. For N = 1218 and 4886 we give the lowest energy configuration we have found, 
though we cannot certain this is a global energy minima. Fivefold coordinated charges (points) 
(pentamers) are indicated by large black dots, and sevenfold coordinated charges (septamers) are 
indicated by small red dots. The rest of the charges are sixfold coordinated (hexamers). Colors in 
online version only. 



FIG. 2: Icosadeltahedral configurations. For a) N = 312 {E = 45629.36272) the symmetric 
icosadeltahedral configuration appears to be a global energy minimum, while for N = 1242 b) 
{E = 747107.43183), c) {E = 747106.46027) addition of defects produces a configuration of lower 
energy (though not necessarily the global energy minimum). As in Figure 1 fivefold coordinated 
charges (points) (pentamers) are indicated by large black dots, and sevenfold coordinated charges 
(septamers) are indicated by small red dots. The rest of the charges are sixfold coordinated hex- 
amers. Colors in online version only. 



FIG. 3: Dihedral configurations. For a) N = 137 {E =8499.53449), b) N = 542 {E =139913.69461), 
the symmetric dihedral lattice configuration appears to be the global energy minimum, while for N 
= 2162 c) {E = 2281595.05127), d) {E = 2281587.48735) addition of defects to the lattice produces 
a configuration of lower energy. We give the lowest energy configuration we have found, though we 
cannot certain this is a global energy minima. Fivefold coordinated charges (points) (pentamers) 
are indicated by large black dots, and sevenfold coordinated charges (septamers) are indicated by 
small red dots. The rest of the charges are sixfold coordinated (hexamers). Colors in online version 
only. 



TABLE I: Energy of icosadeltahedral configurations with N < 400. An * indicates a configuration 
of lower energy, tliough not necessarily the global minimum. For each N, here and Table II, 
we tried fifteen runs-random initial configurations followed by a local gradient descent-to find a 
configuration with an energy lower than that of the icosadeltahedral lattice. For N = A2, 92 and 
162 the best known configuration has exactly twelve pentamers (and thus no dislocation defects), 
but does not have icosadeltahedral symmetry. We note that the split larger analogue of the global 
minimum configuration for A'^ = 42 is not the global minimum configuration for N = 162, and 
neither are the split larger analogues of the N = 92 and 162 global minima (data not shown). 



Charg 


cs 


EiKM'gy 


12 


h= 1 k= 


49.1652530580000 


32 


h= 1 k= 1 


412.261274651000 


42 


h= 2 k= 


732.256241038000 
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61720.0233978130 




* w/defects 


61719.3090545160 
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h= 4 k= 3 


65230.0271225660 


392 


h= 5 k= 2 


72546.2583708950 



TABLE II: Energy of icosadeltahedral configurations with N > 400. An * indicates a configuration 
of lower energy, though not necessarily the global minimum. Thus, for N < 792 we also cannot be 
certain that the symmetric icosadeltahedral configurations are the global minima. 
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